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Abstract 

The comparison of J/ij: photoproduction data in the inelastic region with theoret- 
ical predictions based on the NRQCD approach has remained somewhat ambigu- 
ous and controversial, in particular at large values of the inelasticity variable z. 
We study the polar and azimuthal decay angular distribution of J/ij) mesons as 
functions of z and transverse momentum pt. Future measurements of decay an- 
gular distributions at the HERA ep collider will provide a new test of theoretical 
approaches to factorisation between perturbation theory and quarkonium bound- 
state dynamics and shed light on the colour-octet production fraction in various 
regions of z and pt- 
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1 Introduction 

The production of quarkonium in various processes, especially at high-energy colliders 
(for reviews, see [|l],0]), has been the subject of considerable interest during the past few- 
years. New data have been taken at pp, ep and e^e~ colliders, and a wealth of fixed- 
target data also exist. In theory, progress on factorisation between perturbative and the 
quarkonium bound state dynamics has been made. The earlier 'colour-singlet model' 
has been superseded by a consistent and rigorous approach, based on non-relativistic 
QCD (NRQCD) an effective field theory that includes the so-called colour-octet 
mechanisms. On the other hand, the 'colour evaporation' model of the early days of 
quarkonium physics |Q has been revived Despite these developments the range of 
applicability of these approaches to the practical case of charmonium is still subject to 
debate, as is the quantitative verification of factorisation. The problematic aspect is 
that, because the charmonium mass is still not very large with respect to the QCD scale, 
non-factorisable corrections may not be suppressed enough, if the quarkonium is not part 
of an isolated jet, and the expansions in NRQCD may not converge very well. In this 
situation cross checks between various processes, and predictions of observables such as 
quarkonium polarisation and differential cross sections, are crucial in order to assess the 
importance of different quarkonium production mechanisms, as well as the limitations of 
a particular theoretical approach. In this paper we discuss how polar and azimuthal decay 
angular distributions of J/ip, produced by real photons colliding on a proton target in 
the inelastic region pt > I GeV (or, more conventionally, z = p^ ■ Pp/p^ ■ Pp < 0.9), may 
serve this purpose. 

In the NRQCD approach, to which we adhere in this paper, the cross section for 
producing a charmonium state if in a photon-proton collision is written as a sum of 
factorisable terms. 



where n denotes the colour, spin and angular momentum state of an intermediate cc pair 
and fi/^ and fj/p the parton distributions in the photon and the proton, respectively. 
The short-distance cross sections da{i + j ^ cc[n]) can be calculated perturbatively 
in the strong coupling a^. The matrix elements {O!^) = (0|(9^|0) (see for their 
definition) are related to the non-perturbative transition probabilities from the cc state 
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n into the quarkonium H. The magnitude of these probabihties is determined by the 
intrinsic velocity v of the bound state. Thus the above sum is a double expansion in 
and V. 

Within NRQCD the leading term in v to inelastic photoproduction of J/ip comes 
from an intermediate cc pair in a colour-singlet ^5*1 state and coincides with the colour- 
singlet model result. (The notation for the angular momentum configuration is "^^^^Lj 
with 5, L and J denoting spin, orbital and total angular momentum, respectively.) Cross 
sections , polar , and polar and azimuthal 0] decay angular distributions have been 
calculated for the direct-photon contribution, in which case i = 7 and f^/^{x) = 6{1 — x) 
in (|1]). The angular integrated cross section is known to next-to-leading order (NLO) 
in [0. The colour-singlet contribution, including next-to-leading corrections in a^, is 
known to reproduce the unpolarised data adequately. But there is still a considerable 
amount of uncertainty in the normalisation of the theoretical prediction, which arises 
from the value of the charm quark mass and the wave-function at the origin, as well as 
the choice of parton distribution functions and renormalisation/factorisation scale. 

At order ~ 0.05-0.1 relative to the colour- singlet contribution, the J/tp can also 
be produced through intermediate colour-octet ^Si, ^Sq and ^Pj configurations. In the 
inelastic region, they have been considered in |TO|JTT[| for the direct photon contribu- 
tion and in |12[ for resolved photons, in which case the photon participates in the hard 
scattering through its parton content. Colour-octet contributions to the total photopro- 
duction cross section (integrated over all z and pt) are known to next-to-leading order 
r3| . The polarisation of inelastically produced J/ip due to these additional production 



mechanisms, however, has not been calculated so far. 

Because the colour-octet contributions are suppressed as v*, but, in the inelastic 
region, contribute at the same order in as as the colour- singlet contribution, they are of 
interest only if they are enhanced by other factors, either numerical or kinematical. In 
this respect the situation is similar to a certain f ^-correction, which arises already in the 



colour-singlet model |]T4[ and becomes kinematically enhanced at z close to 1. The colour- 
octet production channels are indeed kinematically different from the colour-singlet one, 
because the ^5*0 and ^Pj configurations can be produced through t-channel exchange of a 
gluon already at lowest order in a^. (For the ^5*1 octet this is true for the resolved process.) 
This leads to a significantly enhanced amplitude, in particular in the large-z region. The 
colour-octet contributions to J/ip photoproduction are indeed strongly peaked at large 
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TO|,|Tl|. Such a shape is not supported by the data, which at first sight could lead 



to a rather stringent constraint on the octet matrix elements (Ogfn]), n G {^^oj^Pq} 
and to an inconsistency with the values obtained for these matrix elements from other 
processes. However, the peaked shape of the ^-distribution is derived neglecting the 
energy transfer in the non-perturbative transition cc[n\ ^ J /ip + X. In reality the peak 
may be considerably smeared as a consequence of resumming kinematically enhanced 
higher-order corrections in f ^ and no constraint or inconsistency can be derived from the 
endpoint behaviour of the 2;-distribution at present. As a consequence, the role of octet 
contributions to the direct process remains unclear. The resolved photon contribution, 
on the other hand, could be entirely colour-octet dominated p!^JT6[| . The ^-distribution 
should then begin to rise again at small z, if the colour-octet matrix elements are as 
large as suggested by NRQCD velocity scaling rules fT^JH and fits to hadroproduction 
data. 

Our motivation for considering the decay angular distributions, including all direct 
and resolved production mechanisms, is to provide another observable that can clarify the 
relative importance of colour-octet production in photoproduction in different kinematic 
regions. Many of the above-mentioned uncertainties and difficulties do not affect the 
polarisation yield. For example, the resummation that is necessary in the endpoint 
region may lead to a significant redistribution of da/dz in z, but affects the normalised 
decay angular distributions to a lesser degree, if they do not have a strong z-dependence 
in the region affected by the smearing. We find that some angular coefficients, especially 
those for the azimuthal angle dependence, take essentially different values in the colour- 
singlet and colour-octet processes. A measurement of decay angular distributions would 
therefore provide information on the relevance of colour-octet production, also at z close 
to 1, which is largely independent of normalisation uncertainties. 

The polarisation of the J / ip can be determined by measuring the angular distribution 
of the leptonic decay J /ip l^l^ . To date, experimental measurements of J /ip polarisa- 
tion exist only for diffractive (elastic and proton dissociation) photoproduction |T8|,|T9[, to 
which the inclusive formalism of NRQCD does not apply, and for fixed-target hadropro- 
duction ^0|]. The latter can be compared with predictions obtained in the colour- singlet 



model |21] and NRQCD P2| , p3[ for p^-integrated cross sections. As discussed in the 
experimental finding of no polarisation is only marginally consistent with the NRQCD 
prediction. Photoproduction offers another opportunity to learn about whether the J / ip 
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polarisation carries information on the spin of the heavy quark pair produced at short 
distances, which is expected in theoretical approaches in which spin symmetry is at work. 
With the expected increase in luminosity at the HERA ep collider, polarisation in pho- 
toproduction of J/ ip at different values of z and pt could provide an attractive diagnostic 
tool in addition to the widely discussed polarisation measurement in pp collisions at the 
Fermilab Tevatron [p^-p7|. 

The paper is organised as follows: Section ^ discusses the production mechanisms 
and calculational details regarding decay angular distributions. In Section ^ we pause 
for theoretical considerations that influence our choice of cuts and other parameters in 
the analysis. Section ^ presents results and their discussion, followed by a summary 
in Section ^. Appendix ^ contains the covariant definitions of coordinate systems and 
polarisation vectors and Appendix ^ summarises the density matrices for all subprocesses 
considered in the paper. 



2 Production mechanisms and cross sections 

We assume that the J /if) transverse momentum pt > 1 GeV, in order to suppress the 
diffractive contribution and higher-twist corrections in general. Away from pt = (or 
z = 1) the leading-twist hard subprocesses contributing to inelastic J/ip production can 
be classified as follows: 

1. Direct photon mechanisms. At leading order in the strong coupling constant, 
O(aa^), these are 



1 + 9 
7 + g/g 



cc 



cc 



3q(1) 3q(8) 1c;,(8) 3p(8) 



+ 9, 



3q{8) lq{8) 3p(8) 



(2) 
(3) 



where the initial-state parton originates from the target proton. 
2. Resolved photon mechanisms. At leading order, 0{al), the subprocesses are 



9 + 9 



cc 



3qW 3q(8) lq{8) 3p(8) 



+ 9, 



q + q 



cc 



3q{8) lq(8) 3p(8) 



+ 9, 
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where one of the initial-state partons originates from the photon and the other 
originates from the proton. 



The direct-photon mechanisms dominate in the region z ^ 0.2, whereas resolved-photon 
mechanisms become important in the region z <, 0.2. (These numbers depend on the 
values of the colour-octet matrix elements, as well as on the pt-cut.) At HERA energies, 
photon-quark fusion can contribute about 10%-15% to the cross section at large z. 
Quark-gluon fusion constitutes about 20%-40% of the resolved cross section ok, z <, 0.2 
and becomes more important than gluon-gluon fusion at larger z. Quark-antiquark 
fusion is always completely negligible. 

The above list includes those colour-octet production channels that are suppressed 
by at most relative to the leading colour singlet production channel. The suppression 
of the octet contributions follows from a multipole expansion of the non-perturbative 
transition cc[n] — > J/il) + X. From a '^Pf^ intermediate state, the physical J/i}) state 
can be reached by a single chromoelectric dipole transition, from a ^S'j^'' state by two 
consecutive electric dipole transitions, and from a ^Sf'' state by a chromomagnetic dipole 
transition. Each electric dipole transition brings a factor f ^, and the magnetic dipole 
transition a factor f^. In addition, the hard production vertex for a P-wave cc state 
is suppressed already by f ^ relative to production in an S-wave state. In photon-gluon 
fusion, the ^S'j^^-amplitude is kinematically identical to the ^S'j^^-amplitude. The ^Sf^- 
channel is therefore insignificant for the direct-photon contribution. 

In resolved photon interactions, on the other hand, the ^5'|^''-channel dominates at 
Pt ^ 5 GeV, because it includes a gluon fragmentation component [Q, in both the gluon- 
gluon and gluon-quark fusion contributions, and therefore falls only as 1/pf at large pt- 
The resolved photon amplitudes are identical to those relevant to J/ip production in 
hadron-hadron collisions ||2^;2^ and at HERA energies the relative importance of the 
various contributions as functions of pt is nearly the same as at Tevatron energies. 

The direct-photon mechanisms above all decrease at least as 1/pf at large pt, with 
the exception of 7 + g — > ccpS*!^^] + q. Fragmentation contributions in photon-gluon 
fusion exist only at the next order in a^. They exceed the leading-order contributions at 
Pt ^ 10 GeV ||30| , p!6[| . We therefore conclude that our list includes all important leading- 
twist production mechanisms for all z and as long as pt ^ 10 GeV. 

We expect that higher-twist corrections due to multiple interactions with the proton 
or photon remnant would be suppressed as a power of A^/ (Q^ + p^), where A ^ 1 GeV 
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is a typical QCD scale and Q is one of the scales involved in the bound state dynamics, 
Q ^ mc,mcV, or rricv'^. Since mv'^ ~ A for charmonium and bottomonium, one may 
expect large higher- twist corrections at small pt, when the heavy quark-antiquark pair 
moves parallel with a remnant jet and remains in its hadronization region over a time 
1/A in the quarkonium rest frame.0 

The differential cross section for J/ip production and its subsequent leptonic decay 
J/ip —>■ l^l^ through any of the resolved-photon subprocesses can be written as 

Bii dildzapt J lons^ z{Xi — z) 



rmn 



X — 
Stt 



Pii + Poo + (Pn ~ Poo) cos^ 6' + V2 Re(p*/o) sin 26 cos + p'i_i sin^ cos : 



where Bu is the J /ip l^l branching ratio, s = {p^ + PpY, s = X1X2S and the parton 
distribution of the proton is evaluated at 

xip2 + M2(xi -2;) 

X2 - J r [ii 

SZ[Xi — Z) 

with factorisation scale fip. Here and in the following we use M = 2mc- The variables 
z and Pt are subject to the restriction 

{l-z)isz-M')>p', (9) 

and 

z{sz-M'^) 
SZ — pf — 

The angles 6 and (j) refer to the polar and azimuthal angle of the in the J/ip decay 
with respect to a coordinate system defined in the J/ ip rest frame. (See Appendix |^ for 
details on their definition.) Finally 

p'ly ^ A[zj ^ J/^P{X) + X] A*[zj ^ J/^(A') + X] (11) 



^Some aspects of higher-twist corrections to 7 + g ^ ccpS'J^^] + g have been considered in with 
the surprising conclusion that the higher- twist correction is A^/ (4m^), even at very large pt, rather than 
A^/pj. The term that does not scale as A^/p^ enters in the combination ei(z) + 4di{z), where ei and 
di are certain twist-4 multi-parton correlation functions defined in [ pT| . However, in the approximation 
considered in one finds ei = Adi. If there existed a sign inconsistency in the A^/(4m^)-tcrm 
would disappear and the result conform to our intuition. 
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are density matrix elements for J/ip production, where a summation (average) over the 
spins of X is understood. The kinematical relations for the direct-photon process 
follow from setting i = 7 and f^/^{xi, fip) = ^(1 — xi). 

The polarisation analysis in NRQCD [p5|j23|j3^ is based on the symmetries of the 
NRQCD Lagrangian, of which spin and rotational symmetry are crucial. In electric 
dipole transitions the heavy quark spins remain intact, so that the J/ip spin orientation 
will be the same as the perturbatively calculable orientation of the total cc spin Sq + 
in the intermediate state. The ^S^^^ intermediate state is rotationally invariant and leads 
to random orientation of the J/ip spin. Technically, we have 



Pa'a' = pL' ?S?] + pL' P^f] + pL' i'sP] + Pix' \{S=l.L = + . . . , (12) 



where p'ly^i [n] refers to production through a cc pair in a state n. The above decomposition 
implies that no interference occurs between the amplitudes for the different terms in the 
sum. The symmetries of NRQCD do not forbid interference of different ^Pj-states. One 



p'l^\{S = 1,L= oc Y: m ^ cc[(lL,; lA)] + X] A*[ij cc[(lL,; lA')] + X] 



J=0,l,2 

where the quantum numbers of the cc pair refer to {LLz, SSz). NRQCD factorisation 
implies that the density matrices can be written as 



where (C^av Wafo..-) ^ NRQCD matrix element with Cartesian indices a,6, and 
K'^'' [n]ab... the corresponding short-distance coefficient. The final step is a tensor decom- 
position of these matrix elements, which, in the case of interest, can be formulated as a 
projection of the cc production amplitude. For J/tp production at the considered order 
in f^, the symmetries of NRQCD are sufficient to reduce all non-perturbative input to 
the four parameters {0'^^'^[n\) with n G {'^S'|"'^\^5'|'^V'S'o^''!^-Po^''} defined as for unpolarised 
J/ip production. 

The calculation then consists of evaluating the density matrix elements for each 
separate term in ( [T^ ) and all partonic subprocesses. We express these matrices as 



finds m 



1^ E pL'?pf\ 



(13) 



(14) 



p'iM = Ann] [e*(A)-e(A')] 



7 



+ M^B'^[n] [pi ■ e*(A) pi ■ e(A')] + M^C'^[n] [p2 ■ e*(A) p2 ■ e(A')] 

+ M^D'^ln] [p, ■ e*(A) p2 ■ e(A') +P2 • e*(A) pi ■ e(A')], (15) 

where e(A) is the J/ip polarisation vector, pi is the momentum of the photon (or the 
parton originating from the photon in resolved contributions) , and p2 is the momentum 
of the parton in the target. The coefficients A, B, C, D are independent of the choice 
of axes in the J/ip rest frame and proportional to a NRQCD matrix element. Their 
analytic expressions are collected in Appendix |B[ 

The decay angular distribution in the J/ip rest frame is often parametrised as 

ocl + A(?/)cos^^ + /i(y)sin2^cos0+ ^4r^sin2^cos20, (16) 



dVtdy 

where y stands for a set of variables and A,/i, are obviously related to (appropriate 
integrals of) the density matrix elements as 

, Pii-Poo A/2Repio 2pi 1 

A = , p = , V = ■ . (17) 

Pii + Poo Pii + Poo Pii + Poo 

Because of the dependence of e(A) on the definition of a coordinate system (see Ap- 
pendix 1^), the parameters A, fj,, v depend on this definition. 



3 Theoretical considerations 

In this Section we discuss some theoretical issues that infiuence our choice of cuts. We 
also motivate the values of NRQCD long-distance matrix elements that we subsequently 
use. 

The NRQCD expansion of the quarkonium production cross section applies to the 
leading-twist contribution of an inclusive production cross section. Leading-twist means 
that the result is accurate up to corrections that scale as some power of A/rrj-c in the 
limit that — > 00. Up to such corrections, NRQCD also applies to the total J/ip 
photopro duct ion cross section. The leading contribution is 0{aas) and purely colour- 
octet [ pil|J53[ ]. It formally contributes only aX z = 1, pt = 0, i.e. in the diffractive 
region. Soft-gluon emission during conversion of the colour-octet cc pair into a. J/ip is 
expected to 'smear' the delta-functions at z = 1 and pt = over a region 6z ~ 0.25, 
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5pt ~ 0.5 GeV |]I5|. One may ask whether the experimentally measured diffractive J/ip 
cross section (with or without proton dissociation) could be considered as part of the 
leading-twist total cross section. Or whether it should be considered as a pure higher- 
twist phenomenon which cannot be regarded as dual (in the sense of parton-hadron 
duality) to the 0{aas) contribution in the inclusive formalism. 

In order for the first possibility to be realised, the soft gluons, which are emitted in the 
transition of the colour-octet cc pair into J /ip+X, would have to recombine into a proton 
or a low-mass diffractive final state. Although it cannot be argued from first principles 
against this possibility, it certainly appears unlikely. It would also hardly be compatible 
with the factorisation assumption of NRQCD that the above colour neutralisation is 
universal, i.e. independent of the rest of the process, again up to higher-twist corrections. 
(Clearly, complete independence is not possible, because some colour exchange between 
J/ip + X and the rest of the process is necessary, if the J/'^ + X state originates from a 
colour-octet cc-pair.) 

The clearest indication that the diffractive contribution should be considered as a 
higher-twist correction, which is not part of a leading-twist calculation of NRQCD, is 
experimental. The HI collaboration has measured [|l^ the polar decay angle distribution 
and the ZEUS collaboration has measured [|T9| the polar and azimuthal decay angular 
distribution. Models of diffractive production based on hard two-gluon p4| , |35| or soft- 
pomeron exchange predict A = 1 (A is defined in (|TBp), in agreement with 
the HERA measurements and earlier fixed-target data |3^. On the other hand, the 
polarisation signature of the leading-twist parton reaction '-yg —>■ cc is identical to the 
signature in the process gg — > cc ||2^. The result is A = if the ^5*0 configuration 
dominates and A = 1/2 if ^Pj dominates. Any linear combination of these values is 
incompatible with the experimental data. Since the diffractive cross section (according 
to the experimental definitions of [l^Jl9[| ) is about as large as the inelastic cross section 
T8| , |38| , we conclude that NRQCD cannot be used to predict the photoproduction cross 



section integrated over all z and pt- 

In order to apply NRQCD we therefore have to cut the elastic region without re- 
stricting the inclusive nature of the process. The HERA collaborations conventionally 
define the inelastic region through the requirement z < 0.9. Let us now argue that it is 
theoretically advantageous to define the inelastic region through a cut in pt. It is obvi- 
ous theoretically, and confirmed experimentally, that the slope of the prdistribution is 
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significantly smaller for inelastic production than for elastic production (with or without 
proton dissociation). A pf-cut at pt > 1 GeV already eliminates most of the diffractive 
contribution as well as higher-twist corrections in general and no further cut on z is 
necessary. In fact, the cross section with an additional cut z < 0.9 cannot be reliably 
predicted in NRQCD. As emphasised in [jl5[, because the NRQCD expansion is singular 
at z = 1, only an average cross section over a sufficiently large region close to 2; = 1 can 
be predicted. The z-distribution itself requires additional non-perturbative information 
in the form of so-called shape functions. These shape functions are also required to 
predict the pt-distribution with an additional cut z < 0.9, but not if z is integrated up 
to its kinematic maximum. In the following, we define the inelastic region through the 
cut Pt > 1 GeV. If statistics is not a limitation, it might be preferable to use Pt > 2 GeV 
to further suppress the higher-twist contributions and difficulties in predicting the pt- 
distribution at low pt, because of (perturbative) soft-gluon emission. Note that the 
resummation of higher-order f ^-corrections in NRQCD will also cause some smearing 
in transverse momentum, which we expect to be less important than that caused by 
perturbative soft-gluon emission. 

Because the colour-octet contributions to inelastic J/ip production are strongly en- 
hanced at large z, an immediate consequence of integrating up to Zmax rather than 0.9 is 
that the p^-distribution is now dominated by colour-octet production, as will be discussed 
in more detail below. The suggested importance of the colour-octet mechanisms could be 
further investigated experimentally, if hadronic activity in the vicinity of the J/ ip could 
be detected. If a J/ip is produced through a colour-octet cc pair, we expect it to be 
accompanied by light hadrons more often than if it is produced through a colour- singlet 
pair. 

The cross sections and decay angular distributions depend on four parameters related 
to the probability of the transition cc[n] J/ip + X. The colour-singlet matrix element 
can be related to the J/ip wave-function at the origin. For {O^^^ {^Si)) we use the value 
obtained in from a fit to hadroproduction of J/ip at large pt- Its precise numerical 
value does not infiuence our analysis, because the ^^i-colour-octet channel is important 
only for resolved photons at large transverse momentum. Our predictions do depend 
crucially on {Oi^'^{%)) and {Oi^^(%)), both of which are not very well known. The 
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Scenario 












I 


116 


1.06 


3.0 




0.0 


II 


116 


1.06 


0.0 




1.0 



Table 1: Values of the NRQCD matrix elements in 10 ^ GeV^ taken for the analysis; 



mc = 1.5 GeV. 



following constraints can be obtained from other J/ip production processes, 

(3.90 ± 1.14) ■ 10-2 GeV^ (Tevatron [26]) 



(Oi^'^(%)) + —(Ot^^'CPo)) = 3.0 ■ 10"' GeV^ (fixed-target hadroproduction m 



3.5 

ml 
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{Oi^'^CSo)) + —AOrCP,)) < 2.8 ■ 10-2 GeV^ 



3.6 



[B ^ J/^X) 



where = 1.5 GeV is assumed. For various reasons, all of these determinations should 
probably be considered as uncertain within a factor of 2. The constraint from inclusive 



B decays has been obtained from the leading-order calculation of ||TT|J39|] , setting the 
colour-singlet contribution, whose magnitude is rather uncertain, to zero. (With the 
parameters of [^, we would have obtained 4.4 ■ lO'^ GeV^ instead of 2.8 ■ lO'^ GeV^.) 
Including the colour- singlet contribution would strengthen the inequality considerably, 
but this cannot be justified given the NLO result of In view of these uncertainties 
and given that they do not allow us to constrain separately {O'l^^ (^Sq)) and {O^^^ {^Pq)) 
with confidence, we consider two scenarios in which the constraints are (approximately) 
saturated either by {Oi'^{%)) or {Oi'^^^Po)) alone. 

The values of all parameters are summarised in Table Further constraints could 
be obtained from the prdistribution in photoproduction, if all kinematically allowed z 
are integrated over. However, for the reasons mentioned earlier, no constraint can be 
derived from the endpoint region of the 2;-distribution. 



4 Results 
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4.1 Cross sections 



We begin with difFerential cross sections in order to display the relative magnitude of the 
various contributions, whose different polarisation yield will influence the decay angular 
distributions. 

The J / ip energy distribution is shown in Figure as a function of the scaling variable 
z = ■ Pp/P'y ■ Pp for a typical HERA photon-proton centre-of-mass energy ^/s^ = 
100 GeV and compared with HI and ZEUS data. (Apart from slightly different 
colour-octet matrix elements, the presentation coincides with that of ||12[-) The colour- 
octet contributions exceed the colour-singlet contribution both for large z ^ 0.65 and 
for small z <, 0.25. 

The normalisation of the short-distance cross sections is strongly affected by the 
choice of the charm quark mass, the QCD coupling, the renormalisation/factorisation 
scale fi, and the parton distribution functions. Varying the parameters in the range 
1.35 GeV <mc< 1.65 GeV, mc</i<4mc, and 150 MeV< A^'') < 250 MeV, the normalisa- 
tion of da{ij —>■ cc[n]) is altered by ~ ±50% around the central value at = 1.5 GeV, 
/i = 2mc, and A^^) = 200 MeV.Q Adopting e.g. the MRS(R2) set of parton distribu- 



tions and the corresponding value of as decreases the short-distance cross sections 
by about a factor of 2 as compared to the leading-order GRV parametrisation. How- 
ever, the values of the non-perturbative colour-octet matrix elements as extracted from 
fits to the Tevatron data |^ depend on the choice of rric, as, and the parton dis- 
tribution in approximately the opposite way such as to compensate the change in the 
short- distance cross section. At leading-twist and leading order in a^, the overall nor- 
malisation uncertainty of the colour-octet contributions to J / ip photoproduction is thus 
in the range of only about ±10%, if the short- distance cross sections are multiplied with 
non-perturbative matrix elements that have been extracted from hadropro duct ion data 
using the same set of input parameters. The long-distance factor of the colour- singlet 
cross section {Oi^^ {^Si)) on the other hand can be determined from the leptonic decay 
width and is not very sensitive to the choice of parameters, up to unknown contributions 
from next-to-next-to-leading-order QCD corrections. Consequently, the normalisation 
uncertainty of the short- distance cross section dci{ij —>■ cc[^5'{^'']) is not compensated 
by a change in the long-distance factor and the colour-singlet contribution should be 



•^To study the dependence of the cross section, we use consistently adjusted sets of parton densities 
il|j42|. 
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Figure 1: Colour-singlet (CS) and colour-octet (CO) contributions due to direct and resolved 
photons to the J/V' energy distribution da/dz at the photon-proton centre-of-mass energy 
y's-yp = 100 GeV in comparison with HERA data |1^,38| averaged over the specified range of 
Js^p. The shaded area bounded by the solid lines represents the sum of all contributions accord- 



'7P 

ing to scenarios I and II for the colour-octet matrix elements. The lines corresponding to sep- 
arate colour-octet contributions are plotted for {O'l'''' {^Sq)) = {Ol'''' {^Po)) /ml = 0.008 GeV^. 
The colour-singlet cross section is evaluated in leading order in a<j. Other parameters: 
rric = 1.5 GeV, renormalisation/factorisation scale /i = 2mc, GRV LO proton and photon 



parton distributions |W1[, A 



,(4) 



200 MeV. 
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considered uncertain within a factor of two. Next-to-leading order QCD corrections 
increase the colour-singlet cross section by ~ 20% — 40%, depending in detail on the 
choice of parameters, but do not affect the shape of the J / energy distribution. 

Given the large normalisation uncertainties in particular of the colour- singlet contri- 
bution, no conclusive statement about the size of the colour-octet matrix elements can 
be derived from the J/ip energy distribution in the region z <, 0.8. On the other hand, 
the dramatic increase of the colour-octet cross section at larger z is not supported by the 
data. One should not interpret this discrepancy as a failure of the NRQCD theory itself, 
but rather as an artefact of our leading-order approximation in and f ^ for the colour- 
octet contributions. Close to the boundary of phase space, for z ^ 0.75, the shape of the 
^-distribution cannot be predicted without resumming singular higher-order terms in the 
velocity expansion |T^. This difficulty is exactly analogous to the well-known problem of 
extracting the CKM matrix element \Vub\ from the endpoint region of the lepton energy 
distribution in semileptonic B decay. To constrain the colour-octet contributions from 
the J/ ip z-distribution, the distribution would have to be averaged close to the endpoint 
over a region much larger than f ^ ~ 0.25. 

The low- 2; region is not expected to be sensitive to higher-order terms in the velocity 
expansion. Therefore, if the data could be extended to the low- 2; region, an important 
resolved photon contribution should be visible, if the colour-octet matrix elements are 
not significantly smaller than assumed in Table 

The prdistribution for inelastically produced J/ if) is shown in Figure ^ with a lower 
2;-cut: z > 0.1. As discussed in Section ^ no upper cut in z is necessary or advisable to 
suppress the diffractive contribution, if the transverse momentum is above about 1 GeV. 
With this definition the differential cross section is dominated by colour-octet contri- 
butions, which exceed the colour-singlet contribution by almost an order of magnitude, 
similar to their significance in hadron-hadron collisions at fixed-target energies [^. Ex- 



perimental data from HERA exist only for pt < 3 GeV ||18| , |38| . The data are presented 
with a cut at 2; < 0.9, in which case the differential cross section at = 1 GeV (3 GeV) is 
found to be about a factor of 10 (2) smaller than in Figure 0. The transverse momentum 
distribution at 2; < 0.9 is adequately accounted for by the colour- singlet channel, includ- 
ing next-to-leading-order corrections in [^]. Diagrams with t-channel gluon exchange 
lead to large i^-factors that increase with increasing transverse momentum and harden 
the pt-spectrum of the colour-singlet channel at NLO considerably. We do not expect 
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Figure 2: Colour-singlet (CS) and colour-octet (CO) contributions due to direct and resolved 
photons to the J/ip transverse momentum distribution da/dpt at the photon-proton centre-of- 
mass energy ^s-yp = 100 GeV; z is integrated to its upper kinematic limit. Other specifications 
are as in Figure |l|. 
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a similar strong impact of next-to-leading order QCD corrections on the transverse mo- 
mentum distribution of the colour-octet cross sections. It would be interesting to learn 
whether including all z can lead to stringent constraints on the size of the colour-octet 
matrix elements. However, in order to obtain an accurate theoretical prediction in the 
lower-pt region, ^ 2-3 GeV, perturbative soft-gluon resummation would have to be 
taken into account. We expect that soft-gluon resummation will be more important for 
the colour-octet processes, because there is no Sudakov form factor for radiation off the 
cc pair in the colour-singlet ^5*1 channel, for which there exists a colour dipole moment 
only. 

4.2 Decay angular distributions 

We now turn to the decay angular distributions, which constitute the main result of 
this work. Below we present the z- and Pf-dependence of the polar and azimuthal decay 
angular distribution parameters A, /x, v defined in (|16D, at a typical HERA centre-of-mass 
energy of = 100 GeV. The quasi-real photons at HERA are actually not mono- 
energetic, but have a distribution in energy given approximately by the Weizsacker- 
Williams approximation. However, in general we have found little energy dependence 
in the energy range relevant to HERA (the only exception being the predictions in the 
recoil frame at z <, 0.3) and thus considered a single energy. 

Since the decay angular distribution parameters are normalised, the dependence on 
parameters that affect the absolute normalisation of cross sections, such as the charm 
quark mass, strong coupling, the renormalisation/factorisation scale and parton distri- 
bution, cancels to a large extent and does not constitute a significant uncertainty. 

The parameters A,/i, as function of z are shown in Figures which include 
direct and resolved photon contributions. We computed the decay angular distribution 
parameters in four commonly used frames (recoil or s-channel helicity frame, Gottfried- 
Jackson frame, target frame and Collins-Soper frame) defined in Appendix 0. Each 
plot exhibits the result from the colour-singlet channel alone (dashed line) and the result 
after including the colour-octet contributions. The two solid lines correspond to the two 
scenarios for the colour-octet matrix elements discussed in Section ^. Recall that the 
J/ip is unpolarised, if it originates from a cc pair in a ^Sq state. Thus, in scenario I 
the angular parameters A, fi, v tend to zero in regions where the colour-octet processes 
dominate. 
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angular parameter ^(z) (Vs^p =100 GeV, Pj > 1 GeV) 
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Figure 3: Angular parameter A of the decay angular distribution as a function of z. Direct 
and resolved photons are included. The dashed line is the colour-singlet model prediction. The 
shaded area shows the NRQCD prediction bounded by the choice of parameters according to 
scenarios I and II. Other parameters are as in Figure 

Inspecting Figure ^ we note that in the recoil, target and Collins-Soper frames A 
differs from the colour-singlet prediction only in the endpoint region. The comparison 
looks different in the Gottfried- Jackson frame: for z <^ 0.5 the colour-singlet channel 
yields large and positive values of A, while the colour-octet contributions yield almost 
unpolarised J/ip. The azimuthal parameter /z (Figure ^ turns out to be least interesting. 
We find that in all frames /z is relatively flat and close to zero, for both the colour- singlet 
and colour-octet contributions. The parameter u, on the other hand, is very different in 
the colour-singlet channel and after inclusion of colour-octet contributions, even in the 
intermediate region of z, where the colour-singlet channel dominates. As can be seen 
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angular parameter |i(z) (Vs p = 100 GeV, > 1 GeV) 
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Figure 4: Angular parameter /i of the decay angular distribution as a function of z. Direct 
and resolved photons are included. The dashed line is the colour-singlet model prediction. The 
shaded area shows the NRQCD prediction bounded by the choice of parameters according to 
scenarios I and II. Other parameters are as in Figure 

from Figure ^, this difference is present in all frames and seems to make u the most 
useful parameter to find out about the relative magnitude of colour- singlet and colour- 
octet contributions experimentally. To determine u one could measure the decay angular 
distribution integrated over the polar angle (cf. 



— — cx 1 + -— + -— cos 2(j), 
acpay 6 6 



(18) 



or project on v as follows: 



8 / dVt co?,{2(j)) 



da 
dQdy 



3 /rffi(l-|cos2^)^ 



(19) 
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angular parameter v(z) (Vs p =100 GeV, > 1 GeV) 
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Figure 5: Angular parameter v of the decay angular distribution as a function of z. Direct 
and resolved photons are included. The dashed line is the colour-singlet model prediction. The 
shaded area shows the NRQCD prediction bounded by the choice of parameters according to 
scenarios I and II. Other parameters are as in Figure 



A distinctive signature of colour-octet contributions in the large- 2; region could be of 
interest in connection with the difficulties in predicting the total cross section in the 
endpoint region. However, the endpoint region in Figures is not without problems 
either. The higher-order terms in the velocity expansion that need to be resummed close 
to the endpoint lead to a convolution of the 2;-distribution with certain non-perturbative 
shape functions. These shape functions depend on the production channel (^5*1^^ ^5*0*^ 
and ^Pj^'') but they are the same for all density matrix elements in every given produc- 
tion channel. As a consequence, while the energy distribution itself depends on these 
shape functions, the moments in z of the normalised angular parameters depend only on 
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angular parameter Mp^) (^s p = 100 GeV, z > 0.1) 
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Figure 6: Angular parameter A of the decay angular distribution as a function of pt- Direct 
and resolved photons are included. The dashed line is the colour-singlet model prediction. The 
shaded area shows the NRQCD prediction bounded by the choice of parameters according to 
scenarios I and II. Other parameters are as in Figure 

the difference of the shape functions in the various production channels. Since we do ex- 
pect such differences, especially between the colour- singlet and the colour-octet channels 
(due to the different properties with respect to soft gluon radiation), and since we are 
interested in the ^-distribution rather than its moments, the predictions for the angular 
parameters in the endpoint region still depend on these shape functions. However, this 
dependence is strong only if the angular parameter varies strongly in the endpoint region 
and disappears entirely if its distribution is fiat. 

In Figures we present the same analysis for the pt-distribution. We note that 
z is integrated up to its kinematic maximum. As a consequence the cross section is 
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Figure 7: Angular parameter [i of the decay angular distribution as a function of pf Direct 
and resolved photons are included. The dashed line is the colour-singlet model prediction. The 
shaded area shows the NRQCD prediction bounded by the choice of parameters according to 
scenarios I and II. Other parameters are as in Figure 



colour-octet dominated and the colour-singlet contribution plays no role for the solid 
curves. Since the colour-octet cross section is dominated by the large-z region, the solid 
curves are entirely determined by the polarisation yield of octet mechanisms at large z. 
Contrary to the situation of hadroproduction at the Tevatron, where one expects large 



transverse polarisation due to gluon fragmentation into colour-octet cc pairs pj], the 
photoproduction cross section tends to be unpolarised in the region considered here. 
Therefore, the pt-distributions do not discriminate between the theoretical prediction 
based on NRQCD and that of the colour evaporation model, which always predicts 
unpolarised J j^^. 
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Figure 8: Angular parameter of the decay angular distribution as a function of pt- Direct 
and resolved photons are included. The dashed line is the colour-singlet model prediction. The 
shaded area shows the NRQCD prediction bounded by the choice of parameters according to 
scenarios I and II. Other parameters are as in Figure 



On the other hand, the transverse momentum distribution seems to prove very useful 
to determine the relative magnitude of colour- singlet and colour-octet contributions: If 
the cross section is dominated by the colour- singlet channel, large and positive values 
of the polar angular parameter A (Figure ^ are expected for pt ^ 5 GeV in the recoil, 
Gottfried-Jackson and target frames. A similar distinctive difference between colour- 
singlet and colour-octet channels is visible in the azimuthal parameter u as defined in 
the Collins-Soper frame, see Figure ||. 

The unique transverse polarisation signature of gluon fragmentation could possibly 
be made visible in the resolved photon contribution. If the direct photon contribution 
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is reduced by a cut on the high-z region, the resolved cross section is dominated by 
g — » ccpS'i'^''] at large transverse momentum |]TB[. As for hadroproduction pp| , pB| , we 
expect that pt ^ (5-10) GeV is necessary to suppress sufficiently the other colour-octet 
channels. 



5 Conclusion 

We presented an analysis of the full polar and azimuthal decay angular distributions 
of inelastically photo-produced J/ip based on the NRQCD factorisation approach to 
quarkonium production. The primary motivation of this study is to find observables in 
addition to angular-integrated differential cross sections, which are sensitive to different 
theories and models of quarkonium production. A particular emphasis is on clarifying 
the role of colour-octet contributions suggested by NRQCD and other quarkonium pro- 
duction processes in comparison with the colour- singlet model, which can be considered 
as a successful description of photoproduction as far as present, limited, data on energy 
and transverse momentum distributions is concerned. 

Assuming NRQCD long-distance parameters as constrained by other J /ip production 
processes such as in hadroproduction and B decay, we have found that the azimuthal 
decay angle distribution as a function of z or pt would be extremely useful for discrimi- 
nating between the colour-singlet model and NRQCD and, to a lesser extent, the colour 
evaporation model. We also noted that transverse-momentum distributions integrated 
over all energy fraction z are colour-octet dominated and could give meaningful con- 
straints on colour-octet matrix elements from the angular integrated rate as well as the 
decay angle dependence. 

While this paper was being written, Fleming and Mehen presented a study of 
leptoproduction of J/ip in NRQCD, which is complementary to our photoproduction 
analysis. Contrary to photoproduction, the leading-twist partonic diagrams at 0{aas) 
can be sensibly compared with a total leptoproduction cross section for large enough 
photon virtuality . Fleming and Mehen computed the polar angle distribution due to 
these leading-order mechanisms and also find interesting tests of NRQCD. 

Acknowledgements. We would like to thank J.G. Korner for useful discussions. M.B. 
wishes to thank John Collins for supplying a FORTRAN integration routine. 
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A Polarisation frames 



We collect here the covariant expressions for polarisation vectors in the four commonly 
used frames, following (Note that the metric g^u = diag(— 1, 1, 1, 1) is used there.) 

Let be the photon momentum, Pp the proton momentum, P the momentum of 
quarkonium ip and s = {p-y +Pp)^- We define the auxiliary vectors 

A = py + Pp A^' = A^ 

B = p^-Pp B'' = B^' 

where M is the mass. (We take M = 2mc in the analysis. The proton mass will be 
neglected in the following.) Note that A, B are three-vectors in the ip-iesi frame. We 
then define a coordinate system as follows: 

1. Choose Z^' = a.A^" + /3^E^, with = -1. 

2. Define = a^A^ + /^^.-B^ in the plane spanned by A, i?, orthogonal to Z and 
normalised: X ■ Z = 0, = — 1. 

3. Take Y to complete a right-handed coordinate system in the ^/'-rest frame, 

= j^e'^'^^W^XpZ, (22) 

(g0i23 _ rjj_^g ambiguity in a^, Px left in the second step is fixed by 
requiring Y to point in the direction of x {—pp) in the ip-rest frame, which 
requires azPx — (XxPz > 0. 

The four commonly used polarisation frames are then specified by the choice of Z. In 
the recoil (or s-channel helicity) frame, the Z-axis is defined as the direction of the 
ip three-momentum in the hadronic centre-of-mass frame, that is Z = —A/\A\ in the 
ip-rest frame. In the Gottfried- Jackson frame Z = p-y/\p-y\ and in the target frame 
Z = —Pp/\pp\. In the Collins-Soper iiame the Z-axis bisects the angle between p^ and 
(— Pp), i.e. Z oc p-y/\p-y\ + {—pp)/\pp\. (All three-vectors refer to the ip-rest frame.) We 
then find the covariant expressions for the coordinate axes from the following expressions 
for a^^x, Pz,x- 



A■PP^' 
M2 ' 

B ■ PP^' 
M2 ' 



(20) 
(21) 
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Recoil frame: 



Oiz = 1 — Pz^O 



A-PB-P 

a. 



{{A ■ py - M^s){{A ■ py - {B ■ py - m^s) 



_ ^{A . py - MH 

Px — 



Is {{A ■ py - {B ■ py - M^s) 

Gottfried-Jackson frame: 

M 



OLz = A 



A-P + B-P 

{B ■ py + A-PD-P + M-s 



{A-P + B- P)^s {{A ■ py - {B ■ py - M^s) 

^ _ {A - py + A- PB ■ P - M^s 

{A- P + B ■ P)^s{{A- py - {B ■ py - M^s) 

Target frame: 

Oiz = -Pz = 



a. 



A-P-B-P 

(D ■ py - A- PD- p + 



Px 



{A-P-B- P)^s {{A ■ py - {B ■ py - M^s) 

{A- py - A- PB ■ P - M^s 
{A-P-B- P)^s {{A - py - {B - py - M^s) 

CoUins-Soper frame: 

BP ^ A-P 

OLz^ ; Pz 



Px 



s {{A - py - {B ■ py^) ^s{{A-py-{B-p) 

MA-P 

A-py-{B- py) {{A - py - {b - py - m^s) 

MB-P 

^J{{A - py - {B - py) {{A - py - {b - py - m^s) 
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We note that A-P + B-P = {P^ + M^)/z and A-P-B-P = sz. Finally, the polarisation 
vectors are given by 



v2 



(35) 



B Density matrices 

The density matrices for all processes considered in the paper are given in this Appendix. 
The results given for the resolved photon process apply equally to J/tp production in 
hadron-hadron collisions and have been used in [^. The functions F, a, b, c, d below are 



related to A, B, C, D of (|l^) as A = Fa etc., and we suppressed all sub- and superscripts. 
For the partonic process 1 + 2 — > 3 + 4, the Mandelstam invariants are s = {pi + ^2)^, 



B.l Direct-photon processes 



7 + (7 — > cc 
F 



+ 9- 

im{ATifaalel{Oi''^{^Si)) 



27[(s - M2)(t - M2)(m - M2)]2 
a = -{s^ + si + Pf + M^{2s^ + si + 2i^){s +i)- M^{f + si + 
b = -2{s^ + P) 
c = -2{s^ + u^) 
d = -2s^ 



(36) 

(37) 
(38) 
(39) 
(40) 



7 + (7 — > cc 



^Sf^ + g: a, 6, c, d are the same as (|37D-(^), F is multiplied by 

15 {Oi'^^i'S,)) 



'y + g ^ cc 



1q(8) 
On 



+ 9- 



Fa 



2{ArrraayM^^{%)) 
Mi[{s - M^){i - M'^){u - M2)]2 



(41) 
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3p(8) 



xsu{[u' - M\u - M'')f - 2si{u - My + s^t"} 
b = c = d = 

+ g: 



t^[M{s - M^){t - M^){u - M2)]3 

- M'^si{s' + 9>sH+ 16s^? + + 8s^£^ - 3sf® - 2f ) 
+ MH{As^ + 2^sH + 48s^P + AVsH^ + 18s^t^ - 4st® + P) 

- MH{12f + 6AsH + 99s^F + SOs^P + 33sH^ + 3sP + 3t^) 
+ MH{22s^ + 88sH+ lUs^P + 67 sH^ + 14# + 

- M^H{22s^ + 6d>sH + eis^P + 16sP + i^) 
+ 2M^2s£(6s2 + 13s£+ 5?) - 3M^^si{s + f) 

-2(s + ty(2s^ + s2p + 2t'^) 

+ 2M'^{Qf + 19s^t + 2ls^P + 17s^P + 21s^t^ + 28s^P + 22st* + 6F) 

- 2M^(6s^ + Qs^f - Ibs^P - 24s^P - AsH^ + 15sP + 5f^) 
+ 2M6(2s^ - mH - 45s^P - 31 - 3s£^ + i^) 

+ 8M^s£(3s^ + 7s£+ 2?) - 8M^°s£(s + t) 

2p(2s6 + 6s^t + 5s^P - IsH'^ - Qsi^ - 21^) 

+ 2M2p(-6s^ - 9s^f + 2s^P + 21s2p + 22s£^ + Si^) 

- 2M\2s^ + 6s^£ - lls^P - 20s=^f^ + IsH^ + 27s£^ + llf^) 
+ 2M^(6s^ + 15s^£ - I2s^t^ - m'^P + 8s£^ + Ap) 

- AM^{3s^ + QsH - AfP - 3sP - 2t^) 

+ 4M^°(s3 + 2sH - sP - 2P) - 2M^H{s - i) 

2P{s^ + 2sH - 6s^P - lOsH^ - 7sP - 2P) 

- M^(4s^ + lAsH+3Qs^P + 29s'^P + As^P - 33sH^ - 39>sP - Up) 

+ 2M^(6s^ + 2lsH + 43s^F + 43s'V/^ + l3fP - l^sP - 9P) 

- M^(12s^ + A5sH + 92s^P + 78s^P + QsP - 9p) 

+ 2M^s{2s^ + llsH + 22sP + 9P) - M^H{5s^ + 6si + P) 
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7 + ^7 — ^ cc 



3^(8)" 



+ q: 



9M^su 



b = 4 
c = 8 
d = A, 

where e„ is the electric charge of the light quark q. 



7 + g — > cc 



10(8) 

•^0 



+ q: 



Fa 



9Mt{t - M2)2 
b = c = d = 



7 + g — > cc 



3p(8) 



+ q: 



a = i{i - M'^){f + + 2MH + 2M^) 

b = -8{s'^ + si + M^i) 

c = 8{P - M^) 

d = 4(P - 2M^s - MH) 



(49) 
(50) 
(51) 
(52) 
(53) 



(54) 
(55) 



(56) 

(57) 
(58) 
(59) 
(60) 



B.2 Resolved-photon processes 



g + g ^ cc 



+ g: a, 6, c, d are the same as (^-(|0D, F is multiplied by 



5 a.o 



96 ae? 



(61) 



g + g^cc\^sf^] +g: 



28 



144M3[(s-M2)(t-M2)(M-M2)]2 I V ^ ^ ^ / V ; 

a = (s^ + St + - (2^2 + si + 2p) (s + t) + M^{s'^ + si + P) (63) 

6 = 2(s' + P) (64) 

c = 2(s2 + 'u2) (65) 

d = 2s^ (66) 

g + g^cc ['Si'^]+g: 



Fa = 



48MstM[(s - M2)(t - M2)(m - M2)]2 
X {s2(s - M2)2 + s£h(M2 - 2s) + Pe} 

X 



(s^ - M^s + M^)2 - M(2? + 3iii + 2^2)} (67) 
6 = c = ci = (68) 



5f + 5f ^ cc 



3 p(8) 



^ ^ 5(4.C..)3(03^/-^(3Po)) 

4s2£2ii2[7V4-(^ _ 7V4-2)(£ _ 7V4-2)(^ _ 7V4-2)]3 ^ ^ 

a = sM {s^°f + 55^^^ + 14s^P ^ 26s^t^ + 35s^£^ + 35s^f^ + 26s^f + Us^P 

+ 5s^P + sf ° 

- M2(s10 + lOsH + 2Qs^P + 40s^P + 45s^t^ + 44s^P + 45s^t^ 
+ m^f + 26sH^ + lOsP + P) 

+ M^(5s^ + 40s^i + 84s^? + lOGs^P + lOSsH^ + lOSsH^ + 106s¥ 
+ 84s^r + 40s£^ + 5?^) 

- M^(16s^ + 104s^£+ 215s^? + 291s^l^ + 316s^£^ + 29153^^ + 215s^t*' 

+ 104sr + 16t^) 

+ MS(34s^ + 179s^t + 356s^P + 473s^P + ATSsH"^ + 3565^^ + 179st^ + 34t^) 

- Mi°(44s^ + 193s^£ + 346s^P + 4105^^ + 346s^f^ + 193sf^ + 44f«) 
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- M^^(15s^ + mH + 52s'^P + 43sf^ + 15^^) 

+ M^%3s^ + 6sH + QsP + 3P)} (70) 

4s^2 + 2As^H + GSs^^t^ + UAs^P + IQAsH^ + 176s¥ + 176sH^ 
+ 176s^f + 164s^f^ + 124s^? + 68sH^° + 24sf ^ + At^^ 

- M^(20s^^ + lOAs^H + 250s^P + 397s^F + 481s^£^ + 500s^£^ 

+ 500s^£^ + 481s^r + 397s^P + 250s^P + 104sf ° + 20f ^) 

+ M^(40s^° + IQGfi + 278s^P + 285s^F + 206s''t^ + U6sH^ + 2Q6sH^ 

+ 285s^f + 278sH^ + 166st^ + 40f °) 

+ M^(-AOs^ - 97 sH+ 53s^P + 373s^f^ + 647s^£^ + 647s^£^ + 373s=^f^ + 53s^f 

- 97s^ - 40P) 

+ M^(20s^ - 33s^£ - 368s^? - 751s^? - 920s^£^ - 751s^^ - 368 

- 33sF + 20t*) 

+ M^°(-4s^ + 77s^t + 323s^? + 492s^P + 492s-'^t^ + 323sH^ + 77si^ - 4?) 

- M^^{41sH+ 120s^P + U2sH^ + 120sH^ + 41sP) 

+ M'''{8sH + Ids^P + IGs^f^ + 8si^) (71) 

-As^^P - 20 - AOsH^ - 40 fi^ + 8sH^ + 80s^F + 128s^t* + llGs^P 
+ 68sH^° + 24sf ^ + 4p2 

+ M^{20s^P + 56s^F + 24s^£^ - 147s^£^ - 409s4^^ - 599s^F - 571s^P 

- 370s^P - 148sf - 28f ^) 

+ M\4s^° + 20sH - 16s^P - 48ft^ + IBOsH"^ + Glls^P + 1060s^?^ 
+ 1155s^f + 854s^f* + 394si^ + 84p) 

- M''(20s^ + 88s^£ + 48s^F + 12 + 318s^t^ + 863s^t^ + 1195s^P 
+ lOeis^r + 583sP + 140P) 

+ M^(40s^ + 152sH+94:S^P + 38s^P + 290sH* + 631sH^ 
+ 738s'^P + 513sP + UOP) 

- M^^(40f + 129sH+53s^P + 7s^P + 129sH^ + 264sH^ + 266sP + 84^) 
+ M^^(20s^ + 55s^i+2s'^P - 15s^P + 30sH'^ + 76sP + 28P) 

+ M'^-As'' - llsH + 7s^P + 7sH^ - IVsP - Ap) 
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(72) 



+ rOs^f ° + 24sf ^ + 

+ M^{4s^^ + 22s^H+70s^P + 115s^P + 71 fi^ - 119sH^ - SSls^P 
- 552s^f - 512s^^ - 320s^P - 126s£^° - 24f ^) 

+ M^(-20s^° - lOAsH - 296s^P - 459 ft^ - 352sH^ + 73s^P + 558s^i^ 



+ 744s^f + 574s^£^ + 270sP 



60t 



10\ 



a 

b = 4 
c = 8 
d = 4 



216M3sM(t-M2)2 

52 + ^2 + 2M2£ 



{4(£-M2)2-9sii} 



g + q^ cc 



1 0(8) 
'-'0 



Fa 
b 



g + q ^ cc 



3 p(8) 



F = 



5(47ra,)=^((!)3^/'^(^^o)) 
216M£(f-M2)2 



{5^ + ^^} 



+ M%AOf + 199s^t + 533 + 778s^P + 596s^t^ + 51sH^ - mS^?^ 
- ms'^f - 296sP - 8OP) 

+ M^(-40s^ - 197s^i - 506s^P - 6725^^^ - 460s^f^ - 79s^£^ + 138sH^ 
+ 16AsF + 60^) 

+ M^°(20s^ + 107sH+267s^P + 307 s^P + 185sH^ + 56 fP - 30sf® - 24F) 

+ M^^{-As^ - 31sH - 7As^P - 71sH^ - A6sH^ - lOsP + At^) 
+ M^\4sH + 8s^P + 8sH^ + 4#) 

g + q^ cc p-S'f ^1 + q: 



(73) 

(74) 

(75) 
(76) 
(77) 
(78) 

(79) 
(80) 

(81) 
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b = -sls"^ + si+MH) 

c = 8(P - M^) 

d = 4(P - 2M'^s - MH) 



81MHu{s - M2)2 

b = 4 
c = 4 
d = 



& = c = d = 



3p(8) 



_ 2O(47ra.)3(0,^/-^(3Po)) 
27M3s2(s _ M2)3 

a = -s{s - M^)(P + u^ + 2M^s + 2M^) 

b = 8(P + s£+M2s) 

c = 8(P + s£ - 2MH + M^) 

d = 4(s2 + 2si + 2P + M^s - 2714"^^) 
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